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Abstract. We study the inelastic scattering of a probe particle on a Bose-Einstein condensate confined in 
a double- well potential. We identify prominent signatures of the underlying mean-field phase space in the 
scattering signal and derive an analytical expression for the inelastic scattering cross section. 
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1 Introduction 

Bose-Einstein condensed particles (BEC) trapped in op- 
tical potentials have sparked an ever growing research ac- 
tivity which has led to many beautiful results such as 
the direct observation of the Mott-insulator to superfluid 
phase transition ^ and of Anderson localization of mat- 
ter waves [Ills]. More recently, hybrid quantum systems 
moved into focus, ranging from BECs coupled to microme- 
chanical oscillators [4 , over integrated single-atom detec- 
tors based on field ionization from graphene nanotubes [5] , 
to fluorescence imaging on atom chips [6]. Interestingly, 
even the smallest (non-trivial) lattice, the double-well po- 
tential, bears rich physics like the appearance of macro- 
scopic, quantum self-trapped states \J_,^ and resulting 
complex decay scenarios [9l [TQ[[TT| rr2] , nonlinear Landau- 
Zener crossings [13I I14[ [T5], and also provides insight into 
seemingly unrelated effects, like e.g., its interpretation as a 
bosonic Josephson junction [16i. Recent experiments even 
took advantage of the inter-atomic interaction present in 
such systems to generate squeezed, i.e. entangled, states 
[17] that allow to approach the Heisenberg limit in atom 
interferometry p!8| [T9] 

The wealth of observed physical phenomena stems from 
the many-body nature of the bosonic system. While the 
microscopic description of the double well is given by the 
Bose-Hubbard (BH) Hamiltonian, a basic understanding 
can be gained by the discrete Gross-Pitaevskii equation 
(GPE), that constitutes the mean-field limit of the BH 
Hamiltonian, i.e., it describes the double- well physics for 
large particle numbers. In the GPE, the inter-atomic inter- 
actions lead to nonlinear equations of motion and, albeit 
being integrable [20 , the mean- field dynamics of a BEC in 



a double-well potential can become complex: In general, 
there are two distinct phase-space regions, corresponding 
to the experimentally observed [8 , Rabi-like oscillations 
and self-trapping of the condensate, respectively. 

In our present contribution, we investigate inelastic 
matter-wave scattering of a probe particle on a BEC con- 
fined in a double- well potential, implement able, e.g., as a 
hybrid quantum system on an atom chip. We employ the 
quantum mechanical treatment of the scattering process 
that we introduced in Ref. [21] to study scattering on a 
chaotic three-site trap. The main question we will address, 
is whether and how the underlying phase-space structure 
is manifest in the quantum scattering signal. While in 
Ref. [21 the chaoticity of the target required a statisti- 
cal analysis, we anticipate that due to its integrability, 
the double- well dynamics leaves immediate fingerprints in 
the scattering quantities. 

This paper is structured as follows: In the next section 
we introduce the target system, a BEC in a double well, 
and discuss its mean-field dynamics and corresponding 



quantum properties. Sees. |2.2| and |2.3| are devoted to the 
scattering setup and its formal description via the scat- 
tering matrix, respectively. Our results are presented in 
Sec. [3 



2 Model 

2.1 Scattering target 

The scattering target is defined by N ultra-cold bosons in 
a double-well potential, described by the Bose-Hubbard 
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(BH) Hamiltonian [Hill]: 



Hbh = 



j,(t) 






u 



hjihj -l)-fc(5{62 + &2^i). 



(1) 



Here, h^ are the bosonic annihilation (creation) opera- 
tors, and rij = V-hj is the number counting operator at 
site j{j/7 and k parameterize the on-site interaction and 
the tunneling strength, respectively. Experimentally, both 
parameters can be independently controlled via the height 
of the potential barrier and by additional magnetic fields 
[24] that induce Feshbach resonances. Apart from the total 
energy E^ also the total particle number A^ is a constant 
of motion. The set of Fock basis states is thus given by 
{|ni)}, where n\ G {0,1,..., TV} denotes the boson num- 
ber in well one, resulting in a Hilbert space of dimension 
7V + 1. 

It is advantageous to introduce the angular momentum 
operators \^\ 



{b2b\-hbl)/2 
Jz= {kbl-hb\)/2 



Jy = 'iib2b[ - bib^ 



(2) 



which obey the commutation relation of an su{2) Lie- 
algebra. In that representation, Hamiltonian (IT]) can be 
rewritten as Hbh = UJ^ — 2k Jx^ up to a constant term. 
The conservation of the total boson number N corresponds 
to [J^^Hbh] = and the physical interpretation of the 
operators Ji is as follows: Jz measures the particle imbal- 
ance between the wells, while Jy represents the conden- 
sate's momentum, and Jx bears direct information about 
the relative phase of the condensate's fractions in the left 
and right well. 

In the mean-field limit (corresponding to large particle 
numbers N), the dynamics of the condensate is described 
by the discrete Gross-Pitaevskii equation ^25j: The quan- 
tum operators b are then replaced by time-dependent, 
complex amplitudes A) , leading to the Hamiltonian: 



ncp/N 



UN ^ 



\A,t-k\A\A2-.AlA^\ . (3) 



The A^ obey the canonical equations idAj/dt = 9Hgp/9A* 
that yield equations of motion for four real variables. Sub- 
tracting the two constants of motion, the mean-field dy- 
namics effectively acts on a two-dimensional phase space 
and can be mapped on a spin N/2 Bloch sphere, via a 
transformation to classical angular momentum variables 
Ji (obtained from (|3| after replacing the operators hi by 
the mean- field amplitudes Ai). 



^ Numerically, we add a small bias of the order of 10~^ni, 
in order to break the symmetry of Eq.([T]) and thereby avoid 
unstable macroscopic superposition states [23] . 
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Fig. 1. (color online). Left: Mean- field phase space of a BEC in 
a double-well potential given by the discrete Gross-Pitaevskii 
equation (|3}, for i^ = 5. Self-trapping modes (red) coexist with 
Rabi-like oscillations (black). The two regions are separated by 
the separatrix (blue). There are two elliptic fixed points asso- 
ciated with self- trapped solutions. Right: The rescaled energy 
spectra En/N of the BH Hamiltonian for various particle num- 
bers A, and constant control parameter i^ = 5. The separatrix 
energy ^sep is plotted as dash-double-dotted, horizontal, green 
line. 



The dynamics generated by ([3| is governed by the 
control parameter u = UN /2k [26 : For very small inter- 
particle interaction i^ < 1, the BH system is in the Rabi 
regime, where the condensate oscillates between the two 
wells, i.e., the time-averaged population imbalance be- 
tween them is zero (i.e., {Jz)t = 0)- In the other extreme 
of very strong interactions u > A/"^, the system is in the 
self-trapped regime, where the condensate is persistently 
trapped in one of the wells, i.e., the average population 
imbalance between them is non-zero {Jz)t 7^ 0- 

For intermediate values of the control parameter u ^ 1 
(Josephson regime) the phase space is divided by a sep- 
aratrix located at the energy E^ep/N = k{u/2 -\- 1), and 
both the above dynamical modes coexist [20l[26l[27] . In 
the left panel of Fig{l] we show the phase-space structure 
typical for this intermediate regime: For small initial par- 
ticle imbalances (and phase difference Jx ^ 0.5) Rabi-like 
oscillations (black) rule the phase-space. Here, each tra- 
jectory corresponds to a situation where a fraction of the 
condensate oscillates between the two wells, while, on av- 
erage, the bosons are equidistributed, i.e., {Jz)t = 0- In 
this Rabi region of phase space and for a given energy, 
there exists exactly one trajectory. 

In contrast, for sufficiently large initial population im- 
balance, one observes self-trapped trajectories (red). They 
encircle the stable fixed points located in the Northern 
(Southern) hemisphere and correspond to a persistent par- 
ticle imbalance Jz{t) > {Jz{t) < 0) for all times t. For 
a given energy in the self-trapped region of phase space, a 
pair of two solutions with inverse particle imbalance ex- 
ists. The blue line denotes the separatrix. 

This phase-space structure has a clear signature in the 
quantum energy spectrum of Hamiltonian ^ (cf. right 
panel of Fig. IT]): Energy eigenstates below the separa- 
trix energy E^ep are non-degenerate, since their mean- field 
counterpart are uniquely defined trajectories in the Rabi 
region. The appearance of two self-trapped trajectories of 
the same energy is quantum mechanically reflected in the 
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Fig. 2. (color online). Scattering setup: The probe particle is 
injected into a waveguide that passes a condensate confined by 
a double- well potential. In the contact region between waveg- 
uide and site one of the double well, the probe particle and 
the BEC can exchange energy. The scattering quantities of the 
particle measured on exit from the waveguide carry detailed 
information on the state of the condensate. 



fact that energy eigenstates above the separatrix energy 
are pairwise degenerate]^ 

In the present paper, we concentrate on the Joseph- 
son regime of intermediate interactions (li = 5), and large 
numbers of bosons in the double- well potential (TV > 30). 
In the following, eigenstates whose mean-field counterpart 
are trajectories that show Rabi-like (self-trapping) behav- 
ior are called eigenstates in the Rabi (self-trapped) region 
of the spectrum. 

We note that by construction (cf. Footnote IT]) none of 
the eigenstates of ^ corresponds to a macroscopic su- 
perposition state. The creation of such states, their in- 
teraction with an electromagnetic field, and the resulting 
decoherence properties have been studied in detail by [29^ . 



2.2 Scattering setup 

We probe the double-well system by means of the scat- 
tering setup shown in Fig. [2] [21 : A probe particle with 
momentum k moves in a waveguide that is placed in the 
proximity of say, site one of the double well. When the par- 
ticle approaches the condensate, it interacts with the latter 
leading to an exchange of energy. On exit from the waveg- 
uide, typical scattering quantities of the probe particle, 
like e.g., its inelastic scattering probability, are measured. 
A possible experimental implementation of such scenario 
is provided by atom-chip setups, which allow for both, 
trapping of the BEC and (magnetic) guiding of the probe 
particle. 



^ More precisely, the eigenstates of ([1| with energy above 
Esep are quasi-degenerate [23 , 28 what is lifted by the small on- 
site potential (see Footnote IT]) . Nevertheless, the self-trapped 
states can be considered as pairwise degenerate, since their 
splitting is very small compared to all other scales in the spec- 
trum of ([T|. 



The waveguide in our scattering scheme is modeled by 
two semi-infinite tight-binding (TB) leads with hopping 
term J and lattice spacing a = 1. These two leads are 
coupled with strength Jq to the central lead site j = 0, 
which is closest to the condensate. Thus, Jo < J deter- 
mines the effective coupling strength between the leads 
and the projectile-target interaction region (i.e. for Jq = 
the latter is completely isolated) and hence controls 
the width of the scattering resonances. Consequently, the 
probe-particle Hamiltonian reads: 



Hr 



TB 






^iS'+i 



Jo 



1,0 



j=-i,o 



^is'+l 



-h.c, (4) 



.(t) 



with Cj the annihilation (creation) operators of the probe 
particle at site j of the TB lead. The particle's energy in 
the momentum eigenstate \km) is 6^ = — 2Jcos(/c^), with 
corresponding velocity Vm = 2Jsin(/c^) [30 . 

The probe-target interaction H^^t is assumed to be 
of similar type (i.e. short range) as the bosonic inter- 
particle interaction in the condensate. Hence, it takes non- 
vanishing values only when the probe particle is located 
at the central lead site {j = 0) which is closest to site one 
of the double well, and it is proportional to the number of 
bosons at this site: 



H]r, 



a • cjco ' 



) Til , 



(5) 



a > is a parameter that controls the strength of the 
probe-target interaction. 



2.3 Scattering matrix 

Given the total Hamiltonian 



Hu 



H^ 



TB ' 



1 



H^ 



BH 



H. 



(6) 



we can now define the scattering matrix of our problem: 
The condensate is initially prepared in an energy eigen- 
state \Em)^ while the probe particle is injected with an 
energy Cm- Hence, asymptotically far from the interaction 
region, the total energy is S = Em + e^jjThe open chan- 
nels (modes) in the leads are then determined by energy 
conservation and the fact that the condensate's final state 
is among the energy eigenstates {En} of the BH Hamil- 
tonian. Hence, the open modes are characterized by the 
kinetic energy Cn = £ — En oi the outgoing probe parti- 
cle. In other words, for a given total energy f , the basis 
of our scattering problem is completely characterized by 
the Bose-Hubbard eigenstates, i.e. {\en)^\En)} = {\En)}- 
The transmission block of the iS- matrix then reads Ii2l] : 



[STm = Vd 



I J 



(1 -j)[£- Hbh] - Oini + ijv 



V^, (T) 



where 7 = (Jq/J)^, and v is the velocity operator. In 
the eigenbasis of the BH Hamiltonian, both Hbh and 



^ In our calculations we rescale the BH spectrum (and the 
interaction operator) to lie within the bandwidth of the lead, 
and thereby avoid evanescent modes. 
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Fig. 3. (color online). Left: Logarithmically color-coded snap- 
shot of the interaction matrix Q, for N = 100 particles and 
u = 5. The insets show magnifications of the Q-matrix struc- 
ture in the Rabi (left) and self-trapped (right) regime, respec- 
tively. Right: The expectation value {En\ni\En) = Qnn of the 
bosonic number operator fii versus the level index n, for identi- 
cal parameters. The error bars denote the corresponding stan- 
dard deviation an = ^y {En\nl\En) - {En\ni\En)'^ . 



V become diagonal matrices i.e. [HBHJnm = En^nm and 
Vnm = VnSnm- In contrast, Qnm = (£^n|^i|^m) IS, in gen- 
eral (/c 7^ 0), a non-diagonal matrix, what yields a non- 
diagonal S'T-matrix. This, in turn, corresponds to inelastic 
scattering and in that sense a controls the degree of in- 
elasticity in the scattering process. 

As known from scattering theory, the imaginary part 
of the S'-matrix denominator determines the width of the 
resonances. Thus, in the present case, the coupling ratio 7 
controls this width: For example, in the limit 7^0, the 
resonances become ^-like, while they grow with increasing 
7. For 7 = 1, Eq. (JTl coincides with the S'-matrix for 
inelastic electronic scattering in a ID geometry derived in 
[3T] . Throughout the paper, we fix 7 = 0.1, corresponding 
to the intermediate regime of overlapping resonances. 



3 Results 

3.1 Interaction matrix Q 

We start our analysis by a direct inspection of the probe- 
target interaction matrix Q shown in the left panel of 
Fig. [3) for A^ = 100 bosons. The value of the control pa- 
rameter is ii = 5, what corresponds to intermediate inter- 
atomic interaction strengths in the BEC (see Fig. fl]) and 
will be fixed throughout the paper. One immediately rec- 
ognizes that Q is a banded matrix that is divided into two 
halves by a "blurred" spot, located in the vicinity of level 
n = 56. The behavior of the matrix elements in either re- 
gion is quite different: Consider, for example, the diagonal 
elements Qnn = (^nl^il^n) that represent the expecta- 
tion value of the boson number in well one, with respect 
to the energy eigenstate \En) (see full dots in the right 
panel of Fig. [3|. Below the energy level n = 56, Qnn takes 
the constant value N/2, what corresponds to equidistribu- 
tion of the bosons among the two wells. Above this energy 
level, the expectation value of fii begins to oscillate, e.g., 
in state n = 57 only 31 particles occupy site one, while 
in state n = 58, the majority of 69 bosons populates site 
one, and so on. 



Besides the average value Qnn^ the corresponding stan- 
dard deviation an = ^y {En\ni\En) — {En\ni\En)'^ bears 
complementary information on the number fluctuations 
and is represented in Fig. [3] by the error bars. One clearly 
sees that dn is largest in the vicinity of level = 56 and takes 
its smallest values for the maximally localized states. 

The observed behavior is a direct signature of the un- 
derlying mean-field dynamics of the BEC: We recall from 
the discussion in Sec. |2.1[ that for A^ = 100 bosons, the 
separatrix energy £^sep is located around energy level £^56 
(see black solid curve of Fig. [l]), i.e., around the same 
level at which the expectation value of hi starts to os- 
cillate. Hence, eigenstates {\En)} with expectation value 
Qnn = N/2 correspond to those mean-field solutions that 
show, on average, vanishing particle imbalance {Jz)t = 0, 
and thus belong to the Rabi region of the spectrum. The 
associated number fluctuations grow with increasing level 
index and take their maximum value close to the separa- 
trix, where the mean-field oscillations in the J^-direction 
of phase space are largest (see left panel of Fig. [l]). On the 
other hand, states with index n > 56 are alternatingly lo- 
calized on either well and thus belong to the self-trapped 
region. Their fluctuations are reduced with increasing in- 
dex n. In the mean- field phase space, increasing n corre- 
sponds to trajectories that are closer and closer to the two 
elliptic fixed points, where the oscillations in any direction 
are zero. 

Note also the checker-board structure in the interac- 
tion matrix Q that is present in the self-trapped region 
and which implies that the interaction operator can solely 
induce transitions between self-trapped states located at 
the same well. In contrast, in the Rabi region, the odd 
off-diagonals are strongly populated, while the even off- 
diagonals vanish. This is due to the fact that the operator 
hi does not commute with the parity operator P, and 
can solely induce transitions between states with different 
parity (see Appendix A). Since energy eigenstates in the 
Rabi region of the spectrum have alternating parity, the 
even off-diagonal elements, which couple eigenstates with 
same parity, are zeroj^ In both regimes, the value of the 
off-diagonal matrix elements decreases exponentially with 
the distance from the diagonal. 



3.2 Participation number of St 

How is the structure of Q, that encodes the properties of 
the BEC, reflected in the scattering signal? In the follow- 
ing, we will focus on the inelastic part of the scattering 
signal that, in contrast to the elastic part, accounts for 
all final configurations of the target, and thus bears much 
more detailed information. Consequently, we analyze the 
off-diagonal elements of the iST-matrix. 

As a first step, we determine the number of final target 
configurations \En) for the BEC initially prepared in an 



Although this argument strictly holds only for symmetric 
double wells, we have checked that a small on-site bias does 
not affect the results. 
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Fig. 4. (color online). Participation number PN of the St- 
matrix versus the scaled level index m/N, for i^ = 5, 7 = 0.1, 
a = 1.0, and various boson numbers N. The PN is aver- 
aged over the entire accessible energy axis 8. The peak around 
m/N ^ 0.56 reflects the separatrix (compare right panel of 
Fig.[l}. 



eigenstate \Em) {m ^ n). To this end, we define the par- 
ticipation number PN of each column of the St -matrix: 



PN(m) 



Y^n^m l[^T]r 



En7^ml[^T]r 



(8) 



Since the diagonal term [S-^]mm is excluded from the sum- 
mation, this quantity is meaningful only if the inelastic 
part of the scattering signal does not vanish, i.e. if a > 0. 
Then, PN reflects the number of outgoing channels that 
participate in the inelastic process, and hence takes values 
between one and N . 

In Fig.|4| we plot the PN for intermediate probe-target 
interaction a = 1 and various boson numbers N ^ versus 
the scaled level index m/N. All curves approximately as- 
sume the value two, apart from the strong peak in the 
PN located around m/N = 0.56. Recalling the above dis- 
cussion, we relate its position to the separatrix energy 
£^sep- Accordingly, for an initial preparation of the BEC in 
the self-trapped or Rabi region (i.e. m./N ^ 0.56), mainly 
two outgoing channels take part in the inelastic scatter- 
ing, while this number is enhanced for initial preparations 
near the separatrix energy £^sep- 

The observed behavior is a consequence of the shape 
of the interaction matrix Q (see Fig. fl]) that can be under- 
stood from a semiclassical argument: In the self-trapped 
and Rabi regions, the mean- field counterpart of the inter- 
action operator, n{t) = |Ai(t)p, oscillates with one unique 
frequency [32], what, upon quantization, yields one tran- 
sition frequency, i.e., two off-diagonals in the Q-matrix. 
In contrast, close to the separatrix, the two distinct dy- 
namical behaviors approach. Therefore, in a small energy 
window around E^ep^ several frequencies appear in the 
dynamics of n(t), what leads to a larger number of off- 
diagonal elements in the interaction matrix Q. 
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Fig. 5. (color online). Inelastic cross section pj^ versus the to- 
tal energy 8, in units of the tunneling strength k. For better 
visibility we chose A/" = 30 particles but otherwise same pa- 
rameters as in Fig. El Some representative channels are labeled 
by their index m. The black (red) curves belong to prepara- 
tions of the BEC in the Rabi (self-trapped) region of the spec- 
trum, while the blue curves correspond to initial preparations 
in eigenstates with energies close to the separatrix energy. 



3.3 Inelastic scattering cross section 

Besides the participation number PN of the S'x-matrix, an 
experimentally easily accessible quantity is the inelastic 
scattering cross section 



p™(f)=2^|[5TU|^ 



(9) 



which denotes the probability to be scattered (forward 
and backward) from an incident channel |e^) to any other 
channel le^) 7^ \^m)i oi"^ equivalently, from a initially pre- 
pared condensate eigenstate \Em) to a final eigenstate 

\En). 

In Fig. pj we plot the inelastic cross section p^{£) as 
a function of the total energy £^ for all m. That is, each 
curve corresponds to an initial preparation of the BEC 
in a different energy eigenstate \Em)' We use the same 
parameters as in Fig. |4] and - for the sake of clarity - 
choose a smaller system with A^ = 30 bosons. As the 
channel number increases from 771 = 1 to 771 = 15, the 
maximum of p^{£) - referred to in the following as the 
resonance position - monotonically shifts to larger energy 
values (solid black lines), until, around channels tti = 16 
and m = 17 (solid blue lines), the behavior changes. For 
channels tti > 18 (red lines), the resonance position oscil- 
lates, i.e., for p\^ it is located dX £ = 0.36, p\^ becomes 
maximal at £ = 0.51, while pf^ takes its maximal value 
at f = 0.38, and so on. 

This "splitting effect" is the most drastic consequence 
of the underlying mean-field dynamics. In Fig. [5J black 
curves correspond to the Rabi region, while red curves re- 
fer to self-trapping. Quite intuitively, one relates the oscil- 
lations of the resonance positions to the alternating mean 
occupation Qnn (see right panel in Fig.p|. In other words, 
for a pair of curves (like e.g. p\^ and Pif ), the one with 
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lower (higher) resonance position corresponds to a self- 
trapped state localized in well one (two). 

This intuition is corroborated by a simplified expres- 
sion for p^{£)^ that we derive from first-order perturba- 
tion theory with respect to the probe-target interaction 
strength a (see Appendix B), and which holds for the in- 
termediate values of probe-target interaction a and cou- 
pling parameter 7 considered here: 



n 
rin.Born 



{£) c 2v-„ 



o?vL 



v2'jj2 



[(1 - -i){E - Em) - aQmm? + I'v^ 

X [{Em\nl\Em) - {Em\ni\Emf] . (10) 

Neglecting for a moment the (weak) energy dependence 
of the velocity v^, Eq. (10) takes a Lorentzian shape. Ac- 



cordingly, the resonance positions are well approximated 
by En + a/{l—j)Qnn- In the Rabi region, they thus expe- 
rience a constant shift due to the approximately constant 
spacing of the BH energies En (see right panel of Fig. [l]). In 
contrast, self-trapped states are nearly pairwise degener- 
ate and the spacing between adjacent pairs grows linearly. 
In this regime, the oscillation in the occupation number 
Qnn causes the splitting of the resonance positions. 

We note that the only ingredients in (ITo| are the eigen- 
energy Em of the BH system, as well as the expectation 
value Qmm and variance <j^ of the number counting op- 
erator fii in the corresponding eigenstate \Em)' In com- 
parison with the standard deviation an shown in Fig. [3| 
expression ( 10 ) reproduces quite well the overall trend of 
the inelastic cross section pl^{S): In the Rabi as well as in 
the self-trapped region, p^^iS) assumes its largest (small- 
est) values for initial target energies close to (far from) the 
separatrix energy E^ep- 



4 Conclusion 

The properties of a BEC trapped in a double- well poten- 
tial were analyzed via the inelastic, quantum mechanical 
scattering of a probe particle. Traces of the underlying 
mean-field phase space, like the appearance of self-trapped 
solutions, were unambiguously identified in various, exper- 
imentally accessible quantities. Based only on the expec- 
tation value and variance of the probe-target interaction 
operator, an analytical expression was derived that eluci- 
dates the main observations. Finally, the proposed scatter- 
ing setup represents a non- destructive measurement of the 
condensate, what is in contrast to standard techniques, 
like time-of-flight imaging, that necessarily result in the 
destruction of the BEC. 
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A Parity in the Rabi region of the spectrum 

In the Rabi region and for vanishing potential bias, the en- 
ergy eigenstates \Em) (with Em < E^ep) of the BH Hamil- 



tonian have well-defined parity, i.e. P\Em) = ±|£^m), with 
P the parity operator that interchanges the particles be- 
tween the two wells. We now show that the interaction op- 
erator fii cannot induce transitions between energy eigen- 
states with same parity. 

We first calculate the commutator [P, ni]. in the Fock 
basis |ni. A/" — ni) = |ni), where ni denotes the number 
of bosons in well one. In this basis, P\ni) = \N — rii), and 
we have 



[P,ni]|ni) = (niP-(iV-ni)P)|ni) 
= (2n2-iV)P|ni). 



(11) 



Assuming that the two energy eigenstates \Em) and \En) 
(m / n) have same parity, we calculate the corresponding 
off-diagonal element of hi: 

{En\hi\Em) = {En\PhiP\Em) 

= {En\{2h2 - N)PP + hiPP\Em) 

= {En\h2\Em) 

= {En\N - hi\Em) = -{En\hi\Em) ,(12) 

i.e. {En\hi\Em)=0. 



B Perturbation theory 



In this Appendix, we derive Eq. (10). The starting point 



of this calculation is the Born expansion of the 6'T-matrix 
with respect to the perturbation parameter a: We incor- 
porate the diagonal of the Q-matrix, 

into the unperturbed Hamiltonian Hbh^ and treat solely 
the non-diagonal part of the Q-matrix, 



Q = Q 



(14) 



as the perturbation. Hence, we obtain as first-order Born 
approximation of ^ the following expression: 



5'T,Born = ^ V 

+ ija 



1/2. 



7 



{l-j){£ - Hbh) - OiQ ^ ijv 



,1/2 



(l-7)(^ - Hbh) - olQ ^ i^v 



Q 



,1/2 



{l-j){S - Hbh) - Q^Q + ijv 



(15) 



Note that the first term in (15), termed Sdj is irrelevant 



for the inelastic cross section pin, since it is diagonal. We 
can further use the perturbative result (15), to obtain an 



analytical estimate for the inelastic cross section and the 
transmission probability. To this end, we rewrite (15) as: 



>>T,Born 



= s^ 



D 



inB 



(16) 
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where k = a/7 is the rescaled perturbation parameter and 
B = Sdv~^Qv~^ Sd is the rescaled perturbation opera- 
tor. From Eqs. (|9| and (16) one obtains the inelastic cross 
section: 



Pin, Born — ^^ 



ETD TD* TD TD* 

-Dmn-Dr^jyj^ JDrnm-'^rn- 



= 2K\E^\BB^\Em) - {E^\B\Em?] . (17) 

It is interesting to note that, in leading order of hz, p^ g^^^ 
depends only on the variance of the perturbation operator 
B in the eigenstate \Em) of the Bose-Hubbard Hamilto- 
nian|j With the definition of B and some basic algebra we 
obtain 



n 

r in, Born 



2k^v-^ 



2 2 

7 < 



[(1 - 7)(^ - Em) - aQmm? + l^V^m 



(18) 



E' 



^kV^ 



2 2 

rn 



[{l-^){£-Ek)-aQkk?+l^vl 



'km • 



In order to simplify the above formula, we assess the part 
A in the last equation: We have numerically verified that 
for moderate values of a and 7 < 0.2, each of the A-terms 
is equal or smaller than unity (apart from the v^^ its a 
Lorentzian curve, i.e. we have to assure that v^ ^ !)• 
The inelastic cross section is then given by: 



Pin, Born — ^^m 



2 2 



[(1 - 7)(^ - E^) - aQrnm? + 7'^^ 
[{Em\nl\Em) - {Em\ni\Em?] . (19) 
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